The Darboux transformation and algebraic 
deformations of shape-invariant potentials. 



D Gomez-Ullatef , N Kamran $ and R Milson§ 

f Centre de Recherches Mathematiques, Universite de Montreal, (QC) H3C 3J7 
Canada. 

I Department of Mathematics and Statistics, McGill University, Montreal (QC) H3A 
2K6 Canada. 

§ Department of Mathematics and Statistics, Dalhousie University, Halifax (NS) B3H 
3J5 Canada. 

E-mail: ullateScrm . umontreal . ca, nkamranSmath . mcgill . ca, 
milsonOmathstat . dal . ca 

Abstract. We investigate the backward Darboux transformations (addition of a 
lowest bound state) of shape-invariant potentials on the line, and classify the subclass 
of algebraic deformations, those for which the potential and the bound states are 
simple elementary functions. A countable family, m = 0, 1, 2, . . ., of deformations exists 
for each family of shape-invariant potentials. We prove that the m}^ deformation is 
exactly solvable by polynomials, meaning that it leaves invariant an infinite flag of 
polynomial modules VnT^ C "P^^^ C where is a codimension m subspace 

of (1, z, . . . , z"). In particular, we prove that the first (m — 1) algebraic deformation 
of the shape-invariant class is precisely the class of operators preserving the infinite 
flag of exceptional monomial modules Vn ^ = (1,2;^, . . . ,z"). By construction, these 
algebraically deformed Hamiltonians do not have an sl(2) hidden symmetry algebra 
structure. 



PACS numbers: 03.65.Fd, 03.65. Ge 
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1. Introduction. 

The Darboux transformation [1, p. 210] |i2|, a method based on the factorization of 
second-order operators, is an important technique for the exact solution of the one- 
dimensional Schrodinger equation [Hill]. The transformation is also a key concept in 
supersymmetric quantum mechanics [Sj and the theory of integrable systems lEllZllHl- 
From the point of view of spectral theory, a non-singular Darboux transformation 
can be characterized by the following three possibilities [H]. First, a potential with 
a lowest bound state admits a unique forward Darboux transformation, which deletes 
the ground state. Second, a potential admits a 2-parameter family (the energy and a 
shape parameter) of backward Darboux transformations, each of which adds a lowest 
eigenvalue. Third, there is also a 1-parameter family of isospectral transformations, 
corresponding to two critical values of the shape parameter, which neither add nor 
delete bound states [T^. 

The Darboux transformation relates to exact solutions in several ways [llj. 
Generally, the forward transformation deforms a given exactly solvable potential to a 
new, solvable potential form. However, if a parameterized family of exactly solvable 
potentials is shape- invariant jT2], i.e., closed with respect to the forward Darboux 
transformation, then the forward transformation furnishes an explicit description of 
the spectrum and eigenf unctions, rather than a new exactly-solvable form. Therefore, 
to obtain a solvable deformation of a shape-invariant potential it is necessary to employ 
a backward Darboux transformation. This was first done for specific potentials such as 
the harmonic oscillator IS.^, while the general theory was developed in 

However, the general form of the deformed potential features integrals of 
eigenfunctions of the original Hamiltonian — in contrast to the original potential, which 
is an elementary function, with bound states also described by elementary functions. It 
has been noted fl] that only certain discrete values of the energy and shape parameter 
correspond to an algebraic deformation, one where the potential and the bound states 
remain elementary functions. Such forms are attractive from the modeling standpoint, 
and are also important theoretically, since exact results can be obtained even in critical 
conditions, where numerical techniques break down. 

In the present article we explain the discrete nature of algebraic deformations by 
characterizing such deformations in terms of polynomial modules left invariant by a 
second-order differential operator. The invariant module approach is an alternative, 
inherently algebraic, approach to exact solvability — developed originally to treat quasi- 
exactly solvable Hamiltonians [T5jll6| ll7j . In this approach, one considers a Hamiltonian 
that, in a suitable gauge, preserves an infinite flag of polynomial modules 

PO C Pi C ^2 C . . . C P„ C . . . , Vn = {l,Z,...,Z^). (1) 

Such a Hamiltonian has an upper-triangular action, and it is, therefore, algebraically 
diagonalizable. On the line, there are exactly three potential forms whose Hamiltonian, 
in a suitable gauge, preserves ((T)). These are the classical shape-invariant potential 
families: the harmonic oscillator, the Morse jlH], and the hyperbolic Poschl- Teller ^H] 
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potentials. This invariant module approach provides an alternative explanation for the 
exact-solvability of these shape-invariant potentials. 

The details of this approach can be found in [13 Ej, and some generalizations 
in |20j. The method of invariant polynomial flags has also been applied in quantum 
many-body problems [21]. The question of determining which second-order operators 
preserve a finite-dimensional polynomial module has been previously considered in a 
number of papers, including [^1^ . 

We consider the algebraic deformations of the three shape-invariant potential 
families [THll4p24j . and show that, in each case, the m*^ algebraic deformation produces 
a Hamiltonian that, modulo a gauge transformation and a change of variable, preserves 
an infinite flag of deformed polynomial modules 

' m ^ ' m+1 ^ ' m+2 ^ • • • ^ / „ ^ • • • , {■^) 

where each is a certain codimension m subspace of Vn, i-e. the span of n — m + 1 
polynomials of degree n. 

We study the first deformation (m = 1) in some detail, and show that Vn^ is an 
exceptional monomial module [22] , an invariant vector space spanned by 

l,z'^,z^, . . . (the first power is omitted). 

We also show that higher deformations produce more complicated polynomial modules, 
but we do not analyze these modules here. 

We will discuss monomial modules more thoroughly in a forthcoming publication 
At this point, we would like to mention that exceptional monomial modules also 
arise in the context of N-fold supersymmetry, [21] ■ Our emphasis is somewhat different, 
since we are primarily concerned with the interplay between the backward Darboux 
transformation and the class of the operators preserving an infinite fiag of polynomial 
modules. 

We also note that invariance of (P) is generally achieved by expressing the gauge 
hamiltonian as a quadratic combination of those generators of s[(2), realized as first order 
differential operators, which leave invariant the infinite fiag of polynomial subspaces Vn- 
These operators are called Lie-algebraic and the Hamiltonian is said to have a hidden 
s[(2) symmetry algebra. Lie-algebraic potentials in one dimension have been classified 
in ^7]. However, not all exactly-solvable potentials are generated by a hidden symmetry, 
the Coulomb potential being a notable counterexample [22] . Since the potentials studied 
here preserve ([2]) rather than (^, they lack an sl{2) hidden symmetry algebra structure, 
and thus furnish a further indication that the exactly solvable class of potentials is larger 
than the Lie algebraic one [T ^ ITT) [ IT7j . 

This paper is structured as follows. In the next Section we describe the forward 
and backward Darboux transform. In Section [Hj we discuss exactly solvable operators 
and algebraic deformations. We also exhibit the invariant fiags corresponding to the 
algebraic deformations of the three shape-invariant potential families. In Section 4 we 
consider exactly solvable operators that preserve the exceptional monomial module, and 
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demonstrate that these are precisely the first-fold deformations of the shape-invariant 
Hamiltonians discussed in Sectional 

2. Factorization and the Darboux transformation 

Let U{x), X G (— cxD, oo) be a continuous, real- valued function, and let 

be the corresponding formally self-adjoint differential operator. We fix a formal 
eigenfunction 

r[0] = Ao<^, 

and note that every such corresponds to a differential factorization 
r — Ao = 

where 

a = dx + (Tx, a* = -dx + a^, a = - ln0. (3) 

For this reason we shall refer to as the factorization function, and to Aq as the 
factorization energy. 

Commutation of the factors defines a partner differential operator 

f = aa* = -dxx + U, 

where 

Jj = U + 2axx- 

The operators obey the following intertwining relation: 

ar = fa. (4) 

As a consequence, the first-order operator a relates the eigenfunctions of the two 
operators: given 

r[^] = A^, 

we also have 

f[i^] = Xip, ip = a[ip]. (5) 

To give a rigorous treatment of the Darboux transformation [31211; we assume 
that U{x) is bounded from below, and let H be the unique self-adjoint, semi-bounded 
operator corresponding to r. The partner potential U{x) is continuous if and only if 
is non-vanishing. In this case, the spectrum of U{x) is bounded from below, and we let 
H denote the corresponding self-adjoint, semi-bounded operator. Letting A denote the 
closed operator corresponding to a, we have that a* corresponds to the adjoint A*. We 
therefore obtain the following non-formal factorizations: 

n-X = A* A, n-X = AA\ 
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where the compositions are appropriately restricted. 

In particular, A maps DiH), the domain of H, to Difi.). The spectral properties 
of this transformation are governed by one of the following 3 possibilities jH]. 

(i) Forward transformation: cf) is square integrable, and A defines an isomorphism 
between ker A = 0^ and D(T-C). Thus, a forward transformation exists if and only 
if Ti possesses a discrete spectrum, in which case it is unique. The transformed 
spectrum differs from the spectrum of Ti by the removal of A, the lowest eigenvalue. 
The first-order operator a transforms the n^^ bound state of Ti. to the [n — 1)*^* 
bound state of H. 

(ii) Backward transformation: is square integrable, and A defines an embedding 
of D{7i) into Dijti). The range of the embedding is precisely The spectrum 
of H. differs from that of Ti by the addition of a lowest eigenvalue, namely A, with 
the ground state given by cf)^^. The operator a transforms the rfi^ bound state of 
Ti to the {n + 1)*^* bound state of Ti. 

A 1-parameter family of backward transformations exist for every A strictly smaller 
than the spectrum of 7i. To describe the possibilities, let 0+ and 0_ be the unique 
(up to a multiple) positive solutions of 

r[0±] = A0±, (6) 

with the property that 0± is square integrable near, respectively, ±oo. The desired 
is of the form 

0=S0++t0_, S,t>0. 

(iii) Isospectral transformation: neither nor 0~^ are square integrable. In this case 
A acts as an isomorphism between D{T-C) and D{Ti.). The operator a transforms 
the n^^ bound state of Ti. to the n^^ bound state of Ti.. Two isospectral Darboux 
transformations exist for every A strictly smaller than the spectrum of Ti: one 
corresponding to = 0+, and the other to = 0„. 

3. Algebraic deformations of shape-invariant potentials. 

3.1. Exact solvability 

We will call a Schrodinger operator 

n = -d^. + u, 

exactly solvable by polynomials if Ti. is equivalent, by a change of variable and a gauge 
transformation, to a second-order operator T that preserves an infinite flag of finite- 
dimensional polynomial modules 

Ml C M2 C M3 C . . . , dimA^„ = n. (7) 
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The exact solvability comes about because T is upper-triangular relative to a basis 
adapted to the above flag, and hence possesses an inflnite list of eigenpolynomials. In 
this paper we will make the hypothesis that such an operator is of the form 

T = P{z)d,-, + Q{z)d-,, + R{z), 

where 

P{z) = P2Z^ + PlZ + Po, Po,Pl,P2^^, 

is a polynomial of degree 2 or less, and where Q{z), R{z) are rational functions. 
We transform the second-order eigenvalue equation 

T[f] = Xf 
into a Schrodinger equation 

n[(j)] = X(j) 

by a change of variables 

x = J {-pyK (8) 

and a gauge transformation 

= eV 



where 



P 



z=z(x) 

f\p-\Q-\P'\ (10) 

U = \P" - \Q' - \P-^ (g - \P'){Q - \P') +R (11) 

z=z{x) 

Let 71 denote the range of the function z{x), and let us assume that neither Q{z) nor 
R{z) have singularities for z ^ TZ. We assume that P{z) is not everywhere positive. If 
it were, then we change P —P. With these assumptions we prove the following 

Proposition 1 Exactly one of the following possibilities holds, according to the number 
and multiplicity of the real roots of P{z): 

(i) The are no real roots. Then, TZ = (—00,00), and U{x) is non-singular. 

(a) There is a double root, p. Then, IZ = (p, 00) or (—00, p), and U{x) is non-singular. 

(Hi) There is a unique, simple root, p. Then, TZ = (— oo,p], or [p, 00). Both z{x) and 
U{x) are even functions. The potential is non-singular if and only if 

Q{p) = lp'{p), or Q{p) = lP'{p). (12) 

(iv) There are two distinct roots pi < p2, and IZ = (— oo,pi], or [p2,oo). Then, both 
z{x) and the potential are even functions. The potential is non-singular if and only 
if holds with, p = pi, or with p = p2, respectively. 

(v) There are two distinct roots pi < p2, and TZ = [pi,p2]- Then, both z{x) and the 
potential are periodic functions, and U{x) has a singularity. 
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Proof. Condition (|12j) follows from In order for the potential to be nonsingular, 

the function {Q — ^P'){Q — |P') must vanish at those roots of P{z) that lie in TZ. The 
singularity occurs in case (v), because both Q{z) and P\z) are first-degree polynomials, 
and hence (fT^ cannot hold for both p = pi and p = p2- 

3.2. Shape-invariant potentials on the line. 

Let us begin by showing, that the class of exactly solvable Hamiltonians that preserve 
the infinite flag (0) consists of shape-invariant potentials on the line. Up to a constant, 
the most general form of a second-order operator that preserves is 

T = P{z)d,, + Q{z)d,, (13) 

where P{z) and Q{z) are arbitrary second and first degree polynomials, respectively. 
Since P„ is invariant with respect to affine transformations of the z variable, we can use 
such transformations to modify P{z) and Q{z) without loss of generality. Additionally, 
by rescahng the x variable, and shifting the spectrum, we reduce our analysis to one of 
the following four canonical forms: 





P{z) 




z{x) 


U{x) 


la 


-1 


2z 


X 


x^ 


lb 


-Az 


Az-2 






II 


-z^ 


-(2A + 3)2+ 1 






III 


z{\-z) 


{A-\)z + l-A 


cosh^(|) 


i(i-A2)sech2(|). 



Table 1. Shape-invariant potentials on the line 



By Proposition ^ the Q{z) in cases la, lb, and III is the most general form for 
which U{x) is non-singular. In case II, we use a translation of the x variable to fix the 
form of Q{z). Thus, we have derived the classical shape- invariant potential families: the 
harmonic oscillator (la) (lb), the Morse potential (II), and the hyperbolic Poschl- Teller 
potentials (III). Each of these potential forms is discussed in more detail below. 

Let us now show that these potentials are shape-invariant by construction. A 
parameterized family of potentials is called shape- invariant if the ground state Darboux 
transformation acts by changing the potential parameters, but leaves the form of the 
potential invariant. For the operators in question, the ground state energy is A = 0, and 
hence the corresponding factorization is given by 



n = (-(-p)^9, + i(Q - ip')(-P)-^) [i-pf^d. + i(Q - ip')(-^)"^) 




where z = z{x) is the change of variable defined in (jH)). The commutation of the factors 
produces a Schrodinger operator 7i, which corresponds to the algebraic operator 

f = P{z)d,, + Q{z)d,, 

where 

Q = P'-Q. (14) 
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For this reason, the forward Darboux transformation for these potentials produces 
another potential (fTT|) of the same form, but with potential parameters modified by 

3.3. Algebraic deformations. 

Now, let us isolate the values of the energy and shape parameters for which the 
backward Darboux transformation of the shape-invariant potentials yields an algebraic 
deformation. We will say that a Darboux transformation is an algebraic deformation, 
when the derivative of ln0 in is either a rational function, or a composition of a 
rational function with an exponential. As per (P)()13p. for the shape-invariant potentials 
under discussion, the factorization function is of the form 

where p{z) is a polynomial, z{x) is an elementary function, and f{z) is a hypergeometric 
function. We will say that such an f{z) is of polynomial type if f'{z)/f{z) is a rational 
function. We observe that 

(ln0(x))' = p'(^(x))/(x) + :^||M^'(x). 

Thus, to obtain an algebraic deformation we must demand that f{z) be of polynomial 
type, with f{z) non-vanishing in the range TZ of z{x) (see Proposition^. 

The regular solution of the confluent hypergeometric equation [211 Sec. 6.1], 

zy"{z) + {c-z)y'{z)~ay{z)=0, (15) 

is given by the confluent hypergeometric function 

y{z) = $(a, c, z) = iFi(a, c; z). 

It can be shown |29j Sec. 6.9] that solutions of polynomial type exist if and only if either 
a or c — a is an integer. These solutions, expressed in terms of generalized Laguerre 
polynomials L'^{z), m = 0,l,2,..., are given below: 

yi{z) =$(a, c;z) a = —m, 

oc U-\z) 

y2{z) = 2;^^"^ $(a — c + 1, 2 — c; z) c — a = 1 -f- m, 

oc z'-^L]^^{z) 
yziz) = e^<^{c — a,c] —z) c — a = —m, 

<^ eU^-i(-z) 
y4{z) = z^~'^e^ ^{1 — a,2 — c, —z) a = m, 

oc z^~'^e^Ll^^'^{—z) 

In general, we have yi{z) = ysiz) and y2{,z) = y^^z) [211 Sec. 6.4]. 

The regular solution of the hypergeometric equation Sec. 6.1], 

z{l - z)nz) + (c - (a + 6 + l)z)nz) - ahf{z) = 0, (17) 
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is given by the Gauss hypergeometric function 

f{z) = F{a, b, c; z) = 2^^1(0, b, c; z). 

Solutions of polynomial type exist if and only if (fT7|) is of so-called degenerate type. This 
means that the monodromy group around one of the regular singular points 0, 1, 00 is 
trivial [211 Sec. 2.2]. These solution, expressed in terms of Jacobi polynomials Pm'^\z), 
are shown below: 

fi{z) = F{a, b, 1 + a + b — c; 1 — z) , a = — m, 

ocPt'^\2z~l) a = b-c-m, (3 = c- 

f^(z) = z^~^{l- zY-''-^F{l~a,l-b,l-a-b + c;l- z) a = m + l, 

oc zf^il - zf P!n'^\2z -1) a = c-b-l-m, (3 = 

fz{z) = (1 — zY~"'~^F{c — a,c — b,l — a — b + c]l — z) c — a = — m, 

oc (1 - zYP!n'^\2z - 1), a = -6 - m, /? = c - 1 

fi{z) = z^~'^F{a + 1 — c,b + 1 — c,l + a + b — c]l — z) c — a = m + 1, 

oc z'^Pt''^\2z -1), a = b-l-m, (3 = 1- 

We will need ()16|) to construct algebraic deformations of the Harmonic oscillator and 
the Morse potential, and ([T^ to construct the algebraic deformations of the hyperbolic 
Poschl- Teller potential. 

3.4- The Harmonic oscillator. 

The general solution of 

- + = A0(x), U^^ix) = x\ (19) 

is given by 

A, Ao, Ai) = (Ao $(i - f , i; x^) + x $(| - |, |, x^)) e"^, (20) 

The n}^ bound state is 

, . ^ . J 0ho(a;; 1 + 2n, 1,0), n even, 

I 0ho(a;; 1 + 2n, 0, 1), n odd, 

where Hn{z) denotes the n^^ Hermite polynomial. For A < 1, the nodeless solutions of 
(fTTHl are given by 



0,„(x;A,l,t), \t\<2 (22) 

U 4'' 

with the extreme values of t corresponding to 0±. This follows from the asymptotic 
properties of $ for large x Sec. 6.13.1]. 

Applying (fTH|) with c = 1/2 and a = l/4 — Ao/4, c.f. ((201), it follows that algebraic 
deformations only occur when the factorization energy Aq is an odd integer. We rule 
out yi{z) and y2{z) because we need Aq < 1. We rule out 2/4(2;), because (E^ shows that 
the corresponding eigenfunction always has node. Thus, we are left with factorization 
functions of the form 



^>ho (a;) = e^'^y^ix'^) = 0ho(a;; -1 - 4m, 1, 0) oc e^H2m{ix). 
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In this way we obtain the following algebraic deformations of the harmonic oscillator: 
Ut\x) = - 2 d^^{\ogH2m{ix)) - 2, m = 0, 1, 2, . . . . (23) 



U 

lOr 




-30L 

Figure 1. Deformations ul\^\x) of the harmonic osciUator for m = 0, 1, 2, 3 

The harmonic oscillator is shape invariant, and so the zeroth deformation is again 
a harmonic oscillator, albeit with a spectral shift. The potentials and eigenfunctions 
of the higher deformations are described in [SniE]- The full, two parameter family of 
backward transformations is discussed in p. 

The new spectral minimum is —1 — 4m, and the new ground state is a multiple of 

The rest of the spectrum is unchanged. The Darboux transformation corresponds to 
the operator 

aC^ = <9x - dec log (e^H2m{ix)^ 

_ Q _^ _ 4:imH2m-i{ix) 
H2m{ix) 

Consequently, the higher bound states are 

V'lS °^ [V^hoj-i] = e''^{H2m{ix)y^p'['\x), j = 1, 2, . . . , 

where 

P^r^{x) = 2{j - l)H2m{ix)Hj_2{x) - 2xH2m{ix)Hj_i{x) 

- AimH2m-i{'i'X)Hj_i{x). 

Thus, the even polynomials 

= P2?(a;), z = x'^, k = l,...,n 

together with = 1, span an invariant {n + l)-dimensional submodule of Vn+m- The 
odd polynomials 

rk{z) = x'^p^^l^{x), k = 0,...,n 
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also span an {n + 1) -dimensional invariant submodule of Vn+m- Therefore, algebraic 
deformations of the harmonic oscillator are exactly solvable by polynomials. 

3.5. The Morse potential. 

The Morse potential ^H] has the form 



f/^„(x) = -(A + i)e-^ + ie-2^ (24) 

The function 

k, C+, C_) = J2 C^e^'^--^'"'HTk -A,1t 2k, e"^') (25) 
is, generically, the general solution of the corresponding Schrodinger equation 

-0"(x) + u,ux)(i){x) = -e<p. 

In the singular case where 1 — 2/c is a non-positive integer, the general solution can be 
given as 

</>_(a;; fc, C+, C_) = C+ e'=^-5^""^(-fc - A, 1 - 2k, e"^) (26) 
+ C_ e-'=^-5<="'$(fc - A, 1 + 2k, e"^), 

where \Ef is the irregular solution of the confluent hypergeometric equation. 

There are no bound states if A < 0, and \A\ bound states otherwise, with the rfi^ 
bound state is 

V^_,„(x) oc 0_(x, A - n, 0, 1) oc e("-^)--^-^L2(^-")(e-^), < n < A. (27) 

We will focus on deformations of potentials with bound states only. For A > the 
spectral minimum is —A^, and hence we must have 

A<\k\. (28) 

In the non-singular case we have by |29t Sec. 6.13.1] 

0+ = 0„,„(x;A;,O,l), (29) 
, / , r(l-2A;) Tik-A)\ 

^- = ^„(.;*.i,-f^fi^). (30) 

In the singular cases, when 2A; — 1 or A + A; is 0, 1, 2, . . ., we have 

0_ = 0_(x;fc,l,O). (31) 

The above hold for k > A. For k < —A, the order of C+, C_ is reversed. 

To obtain algebraic deformations we apply (fTB|) with a = —k — A and c = 1 — 2k, 
and consider the four possible factorization functions 

0^(a;) = e'=^-K^l/,(e-^), ^ = 1,2,3,4, 

in turn. 
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Let m be a non-negative integer. For 0i we need k + A = which by (I28|) imphes 
that A < f . Hence, by (E^) (jSH), 

(t)i{x) = (p^^ix; m - A, 1, 0) = (f)-{x), 

and therefore it generates an isospectral deformation. 

For (f>2{x) "we need 1 — k + A = m, which by ()28|) imphes A < ^ ~ ^- Hence, by 
(EH) (Ei, 

02(a;) = 0n.o(a;; 1 + A - m, 0, 1) = 

and therefore it also generates an isospectral deformation. 
For 03 (x) we need k — 1 — A = m, and hence by (j^^j) 

(j)3{x) = </)^o(a;; m + 1 + A, 1, 0). 

By ()30|) this function is nodeless if and only if 
r(l-2/c) _ r(-l -2A- 2m) 

r{-k-A) ~ r{-l-2A-m) ^ °' ' 
which holds for even m, and fails for odd. 
For 04 (x) we need k + A = —m. Hence, 

04(x) = (f)^,{x; -m - A, 0, 1) = 0+(a;) 

and therefore it generates an isospectral deformation. 

It follows that the only algebraic deformations of the Morse potential corresponding 
to backward transformations, correspond to the factorization function 

0M(a;) = 03(a;) 0, g{m+A+l)x+ie-^-2{l+m+A)^_g-x^^ ^ ^ ^ 

The resulting potentials have the form 

U!^\x) = -{A + l)e-^ + ie-2- - 29.. {log L^'^'^^^^\-e~n) ■ (33) 
The Darboux transformation corresponds to the operator 



= 5. -9.(iog0i™)) 



d.-{l + m + A) + \e-^ + 



where 



Applying to the bound states (P7j) . we infer that the bound states of the deformed 
potential are 



,0-A)x-ie— 



e 



where 

pf'\z) = qj{z)rm{z){z + n - m - 2A - 1) - zqj{z) r^iz) + zqj{z) r'^{z) 
q,{z) =Lf^-'\z). 
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Note that the p^p\z) are polynomials of degree j ' + m + 1, and hence, for every n, the 
polynomials 



1, pt'\z), p^r\^). •••,prw, 

span a codimension m invariant subspace of Vm+n+i{.z). We have demonstrated that 
algebraic deformations of the Morse potential are exactly solvable by polynomials. 



U 




100 




X 



200' 



Figure 2. Algebraic deformations t/'™^(a;) of the Morse potential for A = 2.5 and 
m = 0, 1, 2 and 3. 



3. 6. The hyperbolic Poschl- Teller potential. 

The hyperbolic Poschl- Teller potential ^^1, which includes the class of reflectionless 
1-soliton potentials [32]; has the form 

U,,{x) = \{\-A') sech^d). (34) 

The general solution [23 Sec. 2.9] of the corresponding Schrodinger equation 

can be given as 

k, Co, C^) = cosh(f )^^|Co + 1 + A;, -| + i - A;, 1; - sinh^d )) 

+ C, sinh(f ) + I + A;, -I + I - A;, I; - sinh^d ))}, 

where F{a, b, c; z) also denotes the analytic continuation of the hypergeometric function 
to Re (2;) < 0. For A > 1/2 , the potential dSH) has \A - 1] bound states 

V^pt,i(x), 0<j<A-l. 

The even bound states are 

V'pt,2j(a;) oc 0pt(x; 4 - z - i, 1,0) (35) 

oc cosh(|)5"^p/ ^' ^\coshx). 
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The odd ones are 



A ■ 3 



V'pt,2i+i(a;) oc (/)pt(x; I - i - |, 0, 1) (36) 



oc sinh(|) cosh(|)2 P-^' (coshx). 

We focus on deformations of potentials with bound states only, i.e., A > \. The 
spectral minimum is — — A)"^. For \k\ > \ — A, the nodeless solutions of (fT^ are given 
by 

4, (.-kit) |t|<2ni±izlE(i±i±l) (37) 

with the extreme values of t corresponding to 0± .29^ Sec. 2.3.2, Sec. 2.10]. 
To obtain algebraic deformations we apply p8|) with 

a = \ + k-^, h = \-k-^, c=l-A, (38) 

and consider the four possible factorization functions 

0,(x) = cosh(|)i-^/,(cosh2(f)), I = 1,2,3,4. 

We rule out 02 (a^) and 03 (a;) because these are odd functions, and hence have a node. 
The factorization functions of the form 

0i(x) = 0p,(A,|-i-m,l,O) 

oc cosh(f)^"^pi"^'"^\cosh(x)) 
are nodeless for m > A — ^. The factorization functions of the form 
04(a;) = 0p,(A,-|-i-m,l,O) 

oc cosh(f)^+^pi^^'^\cosh(x)) 

are nodeless for all m = 0, 1, 2, . . .. 

Thus, we see that there are two series of algebraic deformations. In order to study 
deformations for all possible m, we focus on the latter series. The resulting potentials 
have the form 

UiT\^) = + |)(^ + I) sech^d) - 2 d,, (logPir^'^^coshx)) . (39) 

The Darboux transformation corresponds to the operator 

«pr^ = <9x - 9^ log 04 (a;) 
Applying Op^^ to the even bound state functions yields 

= i^^"^\x) s^-\w) (40) 



where 



cosh(|) sinh(x) 



{w + l)qmiw) 
w = 2z + I = cosh(a;) 

sf^\w) = {w + l) {p'j{w)q„,{w) - q',^{w)pj{w)} ~ Aq^{w)pj{w), 
qm{w) =P^-''^\w) 
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Figure 3. Algebraic deformations C/p™'' (x) of the liyperbolic Poschl- Teller potential 
with A = A and m = 0, 1, 2 and 3. 



Hence, the deformed operator, conjugated by the gauge factor jj,^"^\x), preserves 
the codimension m submodule of Vm+n{z) spanned by s^^\w), j = 0,...,n. A 
similar result holds for the deformation of the odd bound states (jHBj) . Therefore, 
algebraic deformations of the hyperbolic Poschl- Teller potential are exactly solvable 
by polynomials. 

4. Exceptional monomial modules. 

In this section we characterize the algebraic structure of the first-fold deformations 
(m = 1) described in the preceding section. We will show that this is precisely the class 
of exactly solvable operators that preserves the infinite flag of polynomial modules 

P^'^ c ^ c C . . . C C . . . , (41) 

where 

= span{l, z\ z^..., z-}, vj,'^ = span{l}. (42) 

Will will call such modules exceptional monomial modules. They are exceptional in 
the sense that the family of second order operators that leave them invariant is very 
rich [221123 • To this effect, let us begin by the following 

Proposition 2 A second-order differential operator preserves Vn^ if and only if it is a 
linear combination of the following 7 operators: 



^2 


= z^d,, + 2(1 - n)z'^d, + n(n - l)z' , 


(43) 


^2 


= z^d,, -{n- \)z^d,, 


(44) 


^2 


= z'd,, , 


(45) 


^2 


= zdzz - dz , 


(46) 


rj.i-2) 

^2 


= (9^^ - 2z'^dz , 


(47) 
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=zd,, (48) 
T(i°^ =1. (49) 

A proof of the Propostion is found in [221123 • If the hnear combination contains the 
raising operators Tg"*"^^ and Tg^^'' then the operator will preserve Pn^ but not the whole 
flag (jlH). These cases are called quasi-exactly solvable in the literature [inilinilll! and 
will be analyzed in detail in |2S]- Since we restrict to exactly solvable cases, we shall 
consider only the following linear combination 

T = p2Tf + piT^-^) + poT^-'^ + g^Tf, (50) 

where the additive constant has been neglected. This can be written as 

T = P{z)d,-, + Q{z)d,, (51) 

where 

P{z)=p2Z^ +PIZ + PQ, (52) 

Q{z) = q2 z - pi - 2pQ z'^ , (53) 

are quadratic polynomials whose coefficients p2,Pi,Po and q2 are arbitrary real numbers. 

The exceptional monomial module Vn^ is invariant with respect to scaling of the z 
variable. By also allowing rescaling of the physical variable x, it suffices to consider the 
following canonical cases. 





Piz) 


Q{z) 


z{x) 


la 


{l-z){z + 2^ 


-2A) q2Z + 2A + l-A{l + A)z'^ 


(I + A) cosh(x) - A 


lb 


z{l-z) 


q2Z - 1 


cosh^d) 


Ic 


-{1 + z') 


q2Z + 2z-^ 


sinh X 


Ila 


-{z-iy 


q2Z-2 + 2z-^ 


-(2A + 3)e^ + 1 


lib 


-z" 


q2z 




Ilia 


%(\-z) 


q2Z + 8- 16z-^ 


2x2 + 1 


Illb 


~Az 


q2Z + 1 




IV 


-1 


q2Z + 2z'^ 


X 



Table 2. Second-order operators preserving the exceptional monomial module 



In cases lb, lib, Illb, the operator is of the form shown in (fT!^ . and therefore 
preserves the full Vn and not just the exceptional module Vn^. Thus, these cases 
describe undeformed shape-invariant potentials. Cases Ic and IV correspond to singular 
potentials, and will not be discussed further. 

Proposition n] shows that the non-singular potentials in case la correspond to 

q2 = ±{A+l) 

Both possibilities yield the same potential form, so we take the former. Using (|TT|) (jH^ 
we have 

fccoshx-l , 2y4 + l 



Uix)= -i(A+l)(A + |) + 2 



(coshx-A;)2' 2A + 3' 
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uil\x) + a + 



A\2 

-pt V-^/ ' V4 ' 2' ' 



For similar reasons, for case Ilia potentials we must take q2 = 4. Using (jllll ()23j] 
we obtain 

U{x) = 3 + x^ + 
r(i), 



2x2 + 1 (2,t2 + 1)2' 

For the case Ila, we translate the x variable and do a spectral shift to set q2 = 2^4+3. 
In this way, (HH) §^ yields 

U{x) = (2 + + \e-'^ -{A + |)e-- + ^^^1^, k = 2A + ?>, 

= U^2}{x) + {2 + Af. 

We should note that the above potential form is non-singular only if 2 A + 3 < 0. This is 
unavoidable, in as much as we showed in section ITBl that the odd deformation of Morse 
potentials with bound states produce singular potentials. 

In summary, we have demonstrated that non-singular Hamiltonians that are exactly 
solvable by an infinite flag of exceptional monomial modules and not by the ordinary flag 
are precisely the first-fold algebraic deformations of the non-singular shape-invariant 
potentials. 

Although these new potentials preserve a full flag of polynomial subspaces, and 
therefore are exactly solvable in the sense defined by Turbiner in jTHI, they do not possess 
a hidden s[(2) symmetry algebra structure. This shows that the exactly solvable class 
is wider than the Lie-algebraic one. 
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